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ABSTRACT 

The sudden spin-down event (‘anti-glitch’) observed in AXP IE 2259+586 on 2012 
April 21 was arguably caused by a decay of its internal toroidal magnetic field com¬ 
ponent, which turns a stable prolate configuration into an unstable one. We refine 
previous models of this process by modelling the star’s magnetic field self-consistently 
as a ‘twisted torus’ configuration in non-barotropic equilibrium (which allows us to 
explore a greater range of equilibrium configurations). It is shown that, if the star’s 
magnetic field is purely dipolar, the change in the toroidal field strength required to 
produce an anti-glitch of the observed size can be ~ 10 times larger than previously 
calculated. If the star has a quadrupolar magnetic field component, then an anti-glitch 
of similar magnitude can be produced via a decay of the quadrupole component, in 
addition to a decay of the toroidal component. We show that, if the quadrupole com¬ 
ponent decays, the minimum initial toroidal field strength and the change in toroidal 
field strength needed to produce the observed anti-glitch are lower than in the pure 
dipole twisted torus. In addition, we predict the maximum anti-glitch sizes, assuming 
that they are caused by a change in ellipticity, in four glitching magnetars and discuss 
the implications for energetics of accompanying X-ray bursts. 

Key words: stars: magnetars - stars: magnetic fields - stars: neutron X-rays: 
individual: IE 2259+586 


1 INTRODUCTION 


All neutron stars undergo electromagnetic braking (ILvne & Graham-Smith) (2006j). In addition, some neutron stars are ob¬ 
served to undergo sudden, randomly occurring jumps in ang ular velocity of either sign. Many radio pulsars, as well as 
five magnetarsLj exhibit sudden spin ups, known as ‘glitches’ I Espinoza et al. 2011 : Yu et al. 20131 : Olausen fe Kaspi 2014 : 
Garcia fe Ranea-Sandoval2015l '). One magnetar, the anomalous X-ray pulsar (AXP) IE 2259 +586, has been observ ed to expe- 
rience a sudden spin down , an ‘a nti-glitch’, accompanying an X-ray burst on 2012 April 21 1 Archibald et al. 20131 1. Recently, 
Garcia fe Ranea-Sandoval 1 2015ll proposed that the anti-glitch was caused by a rearrangement of the internal magnetic field 
of the magnetar: long-term decay of the toroidal field component led to a change from a prolate to another less prolate mass 
distribution. In this paper, we present a refinement of the Garcia & Ranea-Sandoval model, by constructing a non-barotropic 
model star in equilibrium, whose density is perturbed by an imposed magnetic field. By considering a non-barotropic model, 
where pressure is not solely a function of density, we relax some constraints on the magnetic field configuration, allowing a 
wider range of astrophysical scenarios to be modelled. 


Magnetars possess some of the strongest known, naturally occurring magnetic fields (Mereghetti, Pons, & Mela tos 2015). 


The external field strength of a neutron star is usually inferred from its spin down. The internal field, however, cannot be 
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meas ured directly. Observations of bursts and giant flares (Ioka 2 00ll : ICorsi fc Qwenll2011fl and precession (iMakishima et al. 


2014 1 in magnetars have been interpreted to indicate that the internal field exceeds the ex ternal field by at least one or¬ 


der of magnitude. Numeri cal simulations favour a ‘twisted torus’ magnetic configuration ( Braithwaite fc Nordlund 20061 : 
Braithwaite fe Spruit 20061 1 . Magnetic fields of these strengths deform the stellar mass distribution away from spherical sym¬ 


metry. 

Recently the magnetic deformation of non-barotropic neutron stars has been calculated. Gravitational wave observations 
(from which we can obtain ellipticity e) can thus be used to infer the internal magnetic struct ure of neutron stars in principle 
( Mastrano et alfeoil ; Mastrano. Laskv. fc Melatosll2013l : Mastrano. Suvor ov, fe Melatodl2 015h. Because the stars are assumed 
to be non-barotropic [e.g., due to entropy or lepton fraction gradients (iReisenegger fe Goldreichl Il992l : iReiseneggerl 1200 ll . 
2009)], a greater range of MHD equilibria can be constructed and analysed, and one can easily construct good analytic 


approximations to the linked poloidal-toroidal twisted torus found in numerical simulations ( Braithwaite fc Nordlund 20061 : 


Braithwaitel 20091 : Akgfin et al. 20131 : Dall’Osso et al. 2015 1. Although magnetars rotate too slowly to be good gravitational- 


wave source candidates, the above models can be used to deduce the internal field strengths of fast-rotating newborn magnetars 
(age 1$ 10 s), which are better candidates for a hypothetical gravitational wave detecti on from distances as far as the Virgo 
cluster i DaH’Osso fc Stella 2007 : DalPOsso. Shore, fe Stella 20091 : Mastrano et al.l[201ll f. 

The AXP IE 2259+586 suffered an overall spin frequency change ~ —5 x 10" 7 Hz over ~ 10 2 d ([Archibald et al. 


2013ll . According to Archibald et al. ( 2013lf . there are two possible explanations for this event: an anti-glitch with /S.v/v = 
—3.1(4) x 10”' followed by a glitch with ISv/v = 2.6(5) x 10”', or a n anti-glitc h with /S.v/v = —6.3(7) x 10” 7 followed by 
another anti-glitch with Aiz/r' = —4.8(5) x 10"'. Bayesian analysis by Hu et ah ( 2014 1 showed that the second explanation, 
a double anti-glitch, fits the data better. Garcia fe Ranea - Sando val (120151 4 interpreted the angular velocity drop in terms of 
a change in e. This interpretation is exciting because it offers a way to infer the internal held strength and structure of AXP 
IE 2259+586 independent of other mea surements o f e, e.g. , f rom f uture gravitational wave observations of similar objects 
at the beginning of their lives. However, Garcia fc Ranea-Sandoval ( 2015h used a uniform-density star as the unmagnetized 
background state, ignored the contribution of the poloidal Held component to e, and assumed a uniformly distributed internal 
Held (rather than some self-consistent, spatially varying equilibrium configuration). We relax these restrictions in three ways: 
we a naly se a non-uniform parabolic density distribution [which is sufficient to approximate an n = 1 polytropic star: see 


Mastrano et al, (20111) and Sec. 2.1 of this paper], a non-barotropic hydromagnetic equilibrium, and a linked poloidal-toroidal 


field structure. We calculate e and explore how much the toroidal field component needs to change to generate the observed 
anti-glitch magnitude. We also explore the possibility that AXP IE 2259+586 has a quadrupolar field component, which can 

contribute to the anti-glitch if it changes. _ 

In Section 2, we recap the non-barotropic deformation calculations in the literature ( Mastrano et al 1l201ll: iMastrano, Suvorov, fc Melatosl 


2015) and apply them to AXP IE 2259+586 for a dipolar twisted torus configuration. We investi gate a ran ge of initi al and 
final states, that give the observed anti-glitch magnitude, comparing our results directly to those of lGarcfa fc Ranea-Sandovall 
( 2015l f . In Section 3, we construct a dipole-plus-quadrupole twisted torus, apply it to AXP IE 2259+586, and repeat the test 
over a range of initial and final states. In Section 4, we calculate the energy released during these changes in field config¬ 
uration and compare them to the observed outburst energy. We keep the dipole poloidal fiel d stren gth B p = 5.9 x 10 13 G 
constant throughout this paper, corresponding to the spin-down-inferred value, as Garcia fc Ranea-Sandoval 1 20151 ) did. In 
other words, like they did, we assume that the frequency change is purely due to a decay of the internal toroidal field (Section 
2) or the quadrupole component (Sectio n 3), and the dipol e po loidal field (which is responsible for overall spin down) remains 
unchanged before and after the glitch ( Garcia fc Ranea-Sandoval 2015l l. In Section 4, we calculate the change in magnetic 
energy accompanying this process and compare it to the observed X-ray burst energy. In Section 5, we predict the maximum 
sizes of possible anti-glitches in other magnetars caused by this mechanism and their plausibility. In Section 6, we summarize 
our results, discuss briefly other possible causes of the anti-glitch, and discuss future work needed to refine the Garcia & 
Ranea-Sandoval anti-glitch model. 


2 DIPOLE TWISTED TORUS 


In section 2.1, we summarize the calculation of lMastran o et ah (120111 '). which relates e to the strength of the toroidal magnetic 
field component in a non-barotropic star with a twisted torus, whose poloidal component is purely dipolar. In Section 2.2, we 
apply the results to the AXP IE 2259+586 anti-glitch. 


2.1 Hydromagnetic equilibrium 


The magnetic field is assumed to be axisymmetric. It i s decomposed into its poloidal and toroidal components and expressed in 
dimensionless spherical polar coor dinates ( r,Q,<j >) as 1 Chandrasekhar 19561 : Mastrano et al. 2011 : Mastrano fc Melatod 2012 : 


Mastrano. Laskv. fc Melatosl 2013l l. 
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AXP IE 2259+586 antiglitch and its internal magnetization 3 


B = B 0 [ripX7a(r,9) x X7(f> + rit/3(a)V(/)\, (1) 

where Bo parametrizes the overall strength of the field, and rjt set the relative strengths of the poloidal and toroidal 
components respectively ( r/p = 1 without loss of generality), a(r, 8) is the poloidal magnetic stream function, and the function 
/3(a) sets the toroidal field component. We consider separable stream functions of the form a(r,8) = f(r)g(9) in this paper; 
for a dipolar configuration, one has a = f(r) sin 2 8. The radial part of the stream function, /(r), is formally arbitrary, as long 
as it results in a field that is axisymmetric, continuous everywhere inside the star (r < 1), continuous with a current-free 
dipolar external field at the surface (r — 1), finite everywhere, and whose current vanishes at r = 1. For analytic simplicity, 
we choose the polynomial (Mastrano et al. 120111 ) 


f(r) = 


35 


r 2 — 


6r 4 3r 6 

T + ~T 


( 2 ) 


The function /3 must be a function of a to ensure that the magnetic fo rce has no azimuthal component, which cannot be 
balan ced in magnetohydrostatic equilibrium given a field of the form (1) dMastrano et al. 2011 L Mastran oj/^askv. fe Melatos 
2013 1. To conform to the numerical simulations of Braithwaite fc Nordlundf l 2006 1 and iBraithwaitel (l2009h . we want the 
toroidal component to be confined to a circumstellar torus near the neutral curve, i.e., the circle where the poloidal field 
component vanishes, located at some radius r = tn and 8 = 7t/2 (with tn = 0.79 for the dipole twisted torus in Sections 2.1 
and 2.2). Hence, in keeping with previous works, we choose 


/3(a) = 


(a — l) 2 for a ^ 1, 


0 


for a < 1. 


( 3 ) 


We treat the magnetic force as a perturbation on a background hydrostatic equilibrium and write the hydromagnetic 
force balance equation as 


go 


-(V x B) x B = V5p + <5pV<I>, 


( 4 ) 


to first order in B 2 /(pop), in the Cowling approximation (5$ = 0), where p is the zeroth-order pressure, p is the z eroth - 
order density, $ is the gravitational potential, and Sp , 8p , 5$ are perturbations of the latter three quantities. Yoshida ( 2013l l 
calcul ated e for a pure d ipole without taking the Cowling approximation. His values of |e| are at most ~ 2 times those found 
bv lMastrano et all ( 20111 ) (depending on the choice of zeroth-order density profile). A thorough, full-perturbation calculation 


without the Cowling approximation is beyond the scope of this paper. 

We do not assume a barotropic star, nor do we solve the Grad-Shafranov equation, so Sp does not have to be a function 
solely of Sp, and therefore the equation of state imposes no restrictions on the field structure. Neutron star mat ter con ¬ 
sists of multiple species which reach a stably stratified, hydromagnetic equilibrium within a few Alfven time-scales llPethick 
19921 : Reisenegger fe Goldreich 1992 : Reisenegger 200lh . This system is not in full chemical equilibrium, however; the rel¬ 
ative abundances of the constituent particles change by weak nuclear interactions [time-scale ~ 10 5 (T/10 8 K)~ 6 yr] and 
diffu sive processes [time-scale ~ 10 9 (B/10 11 T)~ 2 (T/10 8 K)~ 6 yr]. Between the Alfven and the weak nuclear time-scales 
llHovos. Reisenegger. fe Valdivia 2008), the star is in a hydromagnetic equilibrium state, in which the composition is no t 
determined solely by the density or pressure, and density and pressure do not correspond one-to-one ( Mastrano et al.| 20111 ). 

id@ 


Since most magnetars are 1-10 kyr old (as inferred from spin and spin down), the non-barotropic assumption is valid 
For the zeroth-order density, we choose a parabolic density profile 


p{r) = p c (l-r 2 ), (5) 

where p c = l5M,/(8nRl) is the density at origin, M, is t he stellar mass, and R * is the stellar radius. This choice is merely 
for computational simplicity, but, as shown by Mastrano et al. ( 201 ll f . the resulting e agrees within 5 per cent with that from 
a more realistic n = 1 polytrope. 

The stellar cllipticity is given by 


t = 



( 6 ) 


2 Due to the strong inverse dependence of the diffusion and weak nuclear time-scales on temperature, and the rapid cooling of neutron stars 
over ~ 10 2 kyr, chemical equilibrium may never be reached in practice IIYakovlev. Levenfish, &: Shibanov 1999; Pe trovich &: Reisenegger! 
120108 . i201ll ; iGonzalez-Jimenez, Petrovich, fe Reiseneggeij[2015i '). so our non-barotropic assumption may be valid for more than just the 
relatively young magnetars. 
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4 A. Mastrano, A. G. Suvorov, and A. Melatos 

where Io is the moment of inertia of the unperturbed spherical star, the moment-of-inertia tensor is given by 


Ijk = Rl [ d 3 x[p(r) + 5p(r,6)](r 2 5 jk - XjXk), (7) 

Jv 

and the integral is taken over the volume of the star (r ^ 1). The density perturbation 5p is calculated by taking the curl of 
both sides of equation (4) and matching the </>-components: 


dSp 

~d9 


dr 


{V x [(V x B) x B]}*. 


( 8 ) 


poR* d<£> 

For a dipolar twisted torus with a(r, 9) = f(r) sin 2 9, Mastrano et al. ( 201111 and Mastrano fc Melatod ( 2012ll derived the 
following formula relating e to overall field strength, stellar mass, stellar radius, and the relative poloidal and toroidal field 
strengths: 


e = 5.63 x IO -6 


B p 

10 15 G 


R. Vf M , 
10 4 m J \ 1.4 Mq 



(9) 


In equation (9), B p is the surface field strength at the pole (for a pure dipole, one has B p = 2 Bq), A is the ratio of internal 
poloidal field energy to total internal field energy, A = 0 gives a purely toroidal configuration and A = 1 gives a purely poloidal 
configuration. 


2.2 AXP IE 2259+586 anti-glitch 

Armed with equation (9), we repeat the calculation of Garcia fc R.anea-Sandoval ( 2015h . who calculated the change in 


function of the anti-glitch’s frequency change. Section 2 of the latter reference contains the result 

— (10) 

v 3 

where v is the star’s rotation frequency, A v is the change in frequency during the anti-glitch, and ei (ef) is the initial (final) 
ellipticity. Equation (10) is derived by assuming that t he star initially con tains a mostly toroidal field (which induces ei), 
which slowly decays by Ohmic diffusion and Hall drift ( Vigano et al. 20131 4 . Some critical strain is then reached, when the 
crust cracks in a sudden event (i.e., the outburst), and the star settles into a new field configuration with ef > en The change 
in frequency (i.e., the anti-glitch) arises from a change in moment o f i nertia (becaus e angu lar momentum is conserved), which 
is then related to the change in e to give equation (4) in Garcia fc Ranea-Sandoval 1 2015h . reproduce d here as equation (10 ). 

The measu red fre q uency ch a nge a ssociated with the anti-glitch is —5 x 10“' Hz over ~ 10 2 d ( Archibald et al.l 20131 + 
To explain this, Archibald et al. (2013) proposed that the event consisted of (1) an anti-glitch with Av/v = —3.1(4) x 10~' 
followed by a glitch with Av/v = 2.6(5) x 10~ 7 , or (2) an anti-glitch with Av/v = —6.3(7) x 10~ 7 and ano ther anti-glitch 
with Av/v = —4.8(5) x 10 -7 . The second explanation, a double anti-glitch, is statistically favoured (Hu et al. 20141 '). We take 
the measured v alue Av/v = — 6 .3 X 10 - ' [both for definiteness and to facilitate direct comparison with the results presented 
in Section 3 of Garcia fc Ranea-Sandoval ( 20151+ the calculations can easily be repeated for Av/v = —4.8 x 10 7 ] and solve 
for ef in e qu ation (10) for variou s ei. We present the results in Fig s. 1 and 2, to b e com pa red di rectly to Figs. 1 and 2 of 
Garcia & Ranea-Sandoval ( 2015 '). To facilitate comparison with Garcia fc Ranea-Sandoval ( 20151 + we express the toroidal 
field strength as ( Bt ), i.e., Bt averaged over the volume of the torus. We hold B p = 5.3 x 10 13 G constant, corresponding to 
the value inferred from spin down, and allow Bt to change. 

Figure [T] shows th at the toroidal field streng th mus t decrease from its initial value to give an anti-glitch of the observed 
magnitude, as noted by G arcia fc Ranea-Sandov al (2015), with each curve approaching a particular value of ( B)) as ( B{) —► 0. 
This minimum allowed value for (B\) is found by setting ef = e max in equation (10), where the maximum ellipticity e max is 
obtained by setting A = 1 (a purely poloid al configuration) in equation ( 9). However, we require higher initial (B\), as well 
as larger changes in (Bt), than predicted by Garcia fc Ranea -Sandov a| (20151+ For example, i n order for a star with M* = 1.4 


Mq and R * = 10 4 m [solid curves in Fig.[l]in this paper and in Garci a fc Rane a -Sandov a) 


pie^ in oi 

3 120151)1 


to give Au/v — ^6.3 x 10 


the toroidal field must decrease from (B\) = 2 x 10 15 G to (B[) ~ 8.7 x 10 14 G, whereas Garcia fc Ranea-Sandoval ( 2015h 
calculated (B{) « 1.8 x 10 15 G. This is because the poloidal field component tends to deform the star into an oblate shape; 
neglecting it overestimates the prolateness of the star fo r a given (B t ) . This overes t imate of e means that, to obtain a given 
|ei — ef|, we need (Bt) to change more than predicted by Garcia fc Ranea-Sandoval ( 2015 1. 

The above behaviour is displayed from an alternative viewpoint in Fig. [2] We plot — Az//iz as a function of the toroidal 
field strength change A(Bt)/(Bl), where A (Bt) = (Bl) — (B{), for three different values of (B\), with B p = 5.9 x 10 13 G. In 
the same plot, the horizontal bands indicate — Av/v = 6.3(7) x 10~ 7 (upper band, dark grey) and — Aviv = 4.8(5) x 10 - ' 
(lower band, light grey). Comparing the solid curve in Fig. [2] with the dashed curve in Fig. 2 of Garcia fc Ranea-Sandoval 
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<B'f> (G) 


Figure 1. Initial and final toroidal field combination s consisten t with the obser ved anti-glitch in the AXP IE 2259+586, i.e., solutions 
to equation (10) for Az//+ = —6.3 x 10~ 7 [cf. Fig. 1 of lGarcfa fc Ranea-Sandovall (|2015l) ] as a function of initial volume-averaged toroidal 
field strength (J3J), for surface polar field strength B v — 5.9 x 10 13 G (inferred from spin down) and three combinations of stellar masses 
and radii [chosen bv lGarcia &; Ranea-Sandovall (j2015h ]: M* = 1.4Mq, R* = 10 4 m (solid curve); M * = 1.4M©, R * = 1.4 x 10 4 m (dashed 
curve); and M* = 1.8M©, R* = 1.2 x 10 4 m (dotted curve). The thin solid curve is (-BJ) = (-£?/) for reference. 


( 20151) (both representing (B£) = 2 x 1 0 15 G), we see that, to obtain /S.v/v = —6.3 x 10 ', our model requires A (St)/ (B)) « 50 
per cent, whereas Garcia fe Ranea-Sandoval ( 20ld l calculated A {Bt) / {B\) « 5 per cent. Thus, by ignoring the poloidal field 
component, one overestimates |e| in general and underestimates A (Bt) / {B\) for a given | Au/i/\. For example, for 1.9x 10 15 G < 
(Bt) < 4.9 x 10 15 G [correspo nding to the minimum value needed for the anti-glitch and the expected limit of stability; see 
below and Akgfin et al. 1 20131 )] and B p = 5.9 x 10 13 G (inferred from spin d own), we find e/eoR.s ~ 0. 1 5, wh ere e is the 
ellipticity calculated using equ ation (9) and eons is the ellipticity calculated by Garcia^^ane^Sandoya] 1 20151 ). 

The analytic calculation of lAkgiin et al. ( 20131) and the numerical simulation of Braithwaitel 1 20091 ) found that 10~ 3 < A < 
0.8 is required for a stable dipolar twisted torus. For B v = 5. 9 x 10 13 G, th i s corr esponds to 7.7 x 10 13 G < ( B t ) < 4.9 x 10 15 
G. A more stringent limit is set by crust cracking. Horowitz fe Kadau ( 20091 ) conducted large-scale molecular dynamics 
simulations of Coulomb solids to represent a neutron star’s crust. They found that the large density and pressure make the 
crust very strong and rigid, capable of supporting |e| < 4 x 10~ 6 (assuming canonical neutron star mass and radius) before 
cracking suddenly in a collective manner (rather than yielding continuously). If we assume that the initial magnetic field does 
not deform the star sufficiently to crack the crust, then this limit translates to ( B\) < 3.7 x 10 lj G in the case of AXP IE 
2259+586. The values of {Bt) plotted in Figs, [l] and [2] fall within these relatively generous limits. 


3 DIPOLE-PLUS-QUADRUPOLE TWISTED TORUS 


Neut ron star magnetic fields are approximately dipolar at radio emission altitudesjL^n^^^anchester 1988; Chung_&JMelatoi 


2011 ) and in the outer magnetosphere, where high-energy emissions originate ( Romani fe Yadigaroglulll995 : Lvutikov, Otte, fe McCann 
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6 A. Mastrano, A. G. Suvorov, and A. Melatos 



Figure 2. The fractional change in frequency /\u/v as a function of the relative change in the volume-averaged toroidal magnetic 
field strength \A.(Bt)\/ (B\) , where A (Bt) = {B’j) — (of) [cf. Fig. 2 of lGarcia fe Ranea-Sandovall (1 201 5| )1. i.e., solutions to equation (10) 
compared to the observed anti-glitch in the AXP IE 2259+586, for surface polar field strength B p = 5.9 X 10 13 G and three values of (-BJ): 
2 X 10 15 G (solid curve); 3 X 10 15 G (dashed curve); and 4.5 X 10 15 G (dotted curve). The bands show Ar/r for the two possible AXP 
IE 2259+ 586 anti-glitch magnitudes inferred from the 2012 April data [taking the statistically favoured double anti-glitch interpretation 
llHu et al .11*20 1 4 ) and their corresponding error bars: Ai'/v = —6.3(7) X 1 0 7 (the initial anti-glitch, dark grey) and —4.8(5) X 1 0 7 (the 
follow-up anti-glitch, light grey). 


2012). However, some observations can be interpreted as signatures of higher-order multipoles cl ose to the surface iBurnettj^^ehRoj 


2014), e.g. cyclotron resonant scattering line energies of cert ain accretion-powered X-ray pulsars llN ishimural2005|: Privmak. Melat os. fe Laskv 
2014), radio emissions from pulsars beyond the ‘death line’ (Young. Manchester. fc Johnstonll999l: Camilo et all2000t Medin fc Lai 


2010 . the ano malous braking index of some radio pu l sars ( Barsuk ov fc Tsvean 2Q10l) . and the substructures found in some 


pulsar signals i Bonazzola. Mottez. fc Hevvaerti 20151 : Petri 2015 ). While the dipole component of the magnetic field can 
be inferred from the observed spin-down rate, the putative higher-order multipoles contribute small corrections of order 
< (2-7T vRt/c) 2 to the torque and cannot be measur ed directly. 


In this section, we outline the method used by Mastrano. Suvorov, fc Melatos ( 20151 ) to describe a self-consistent, linked 
dipole-plus-quadrupole-plus-toroidal field configuration and to relate the resulting e to (Bl), (£>/), and the amplitude of the 
quadrupolar component. 


3.1 Hydromagnetic equilibrium 

If the poloidal field is a combination of dipole and quadrupole components, then the stream function a(r, 9) can be written 
as a linear combination of the form 


a(r, 9) = /i (r) sin 2 9 + re /2 (r) sin 2 9 cos 6, 


( 11 ) 


where re is a dimensionless parameter controlling the amount of quadrupole field present. In addition to the finiteness and 
continuity conditions discussed in Section 2.1, /i(r) and / 2 (r) must now be chosen such that the Lorentz force and the 
resulting density perturbation are continuous everywhere. This is no t a trivial ta sk and is b est accomplish ed by moving to 
a coordinate system defined by the stream function, as described by Mastrano. Suvorov, fc Melatos! (|2015| ). Suitable fx^ir) 
are then f ound throug h tria l and error. The calculation of e itself is straightforward, using equations (6)-(8). Following 


Mastrano . Suvoro v. fc Melato s ( 20151 ) for simplicity and consistency, we use 


/iM 




r 8 + 



12 . 16 
r + ar , 


( 12 ) 


© ? RAS, MNRAS 000.HHT21 

























































AXP IE 2259+586 antiglitch and its internal magnetization 7 


/ 2 (f) = - (35r 4 - 42r 8 + 15r 12 ) , 
o 


(13) 


where a is a dimensionless free parameter which controls the volume of the torus occupied by the toroidal field component. 
These fi, 2 (r) choices are not unique, but they are among the simplest possible polynomials, with the lowest possible order, 
that guarantee a well-defined field-aligned coordinate transformation. _ 


The stability of dipol e-plus-quadrupole configu r ation s has not been calculated, either analytically [like lAkgiin et al. 


( 2013h ] or numerically [like Braithwaite fc Nordlundi ( 20061 )]. As a consequence, some combinations of a, re, and ( B\'*) may 
be excluded in reality for stability reasons. We defer the numerical calculation (by evolving the field in a time-dependent 
magnetohydrodynamic simulation) to a future paper. 


3.2 AXP IE 2259+586 anti-glitch 


In this section, we repeat the calculation of Gar cia fe Ranea-Sando val ( l2015l l. except that we assume AXP IE 2259+586 
possesses some quadrupolar magnetic component, which partially decayed during the anti-glitch. In other words, we assume 
that the change in e which led to the anti-glitch is due to a change in field geometry, in addition to a decay of the toroidal 
component. We hold Bo = 2.95 x 10 13 G constant, so that the dipolar component’s surface polar strength is 5.9 x 10 13 G (we 
thereby assume that the measured spin down, from which a polar field strength of 5.9 x 10 13 G is inferred, is due solely to 
the dipole component). We also set a = —5 (for the sake of definiteness), and allow re and ( Bt ) to change. 

In Fig. [3] we show (B{) given by equation (10) for Av/v = —6.3 x 10~ 7 [cf. Fig. 1 of Garcia fc Ranea-Sandoval ( 2015l ll 


as a function of (B\), for M* = 1.4 Mq, R * = 10 4 m, Bo = 2.95 x 10 13 G, and three different re/ (final re) values: re/ = 0.6 
(solid curve), re/ = 0.4 (dashed curve), and re/ = 0.1 (dotted curve). We assume re; = 0.8 (initial re), i.e., that the initial 
surface quadrupole magnetic field strength at the pole is 1.6-Bo (surface dipole magnetic field strength at the pole is 2Bo as 
before). This rather extreme value is chosen for the sake of definitenes^]. For the dipole-plus-quadrupole twisted torus, t can 
be written as a function of (Bt) according to 


f Bo V f R- 


Vl0 13 G/ \ !0 4 m/ V 1.4 M 0 


M, 


1 - 


b(B t y 

B 2 0 


(14) 


where the values of a and b are given as functions of re in Table 1 (for re = 0.1, 0.2, 0.4, and 0.8; second and third columns in 
the table). 

Figure [3] tells us that, for a given (B\), (B{) decreases as re/ decreases. This is because, in general, the quadrupole 
component induces more positive (oblate) e: for a given (Bt), a configuration with re = 0.1 is more prolate (i.e., e more 
negative) than one with re = 0.6, for example. Therefore, a decrease in re means a decrease in e, which must be countered by 
a decrease in (Bt) (which increases e) to obtain the (e; — e/) value required to match the observed Avjv through equation 
(10). Compared to the pure dipole case of Section 2, however, we find that (B{) is generally larger for the mixed case than for 
the pure dipole case. For example, with (B\) = 2 x 10 15 G, one finds (B{) = 8.7 x 10 14 G for the pure dipole case (see solid 
curve in Fig. 1), but (B{) = 1.8 x 10 15 G for the re; = 0.8 to re/ = 0.6 transition (solid curve in Fig. 3). Note also that, while 
the general behaviour of the (Bj)-(B/) curves in Fig. 3 are similar to Fig. 1, the minimum values of (B\) are smaller for the 
mixed cases, e.g., (B\) « 3 x 10 14 G for the transitions shown in Fig. 3, compared to 1.8 x 10 15 G for the pure dipole case 
(Fig. 1, solid curve). The decay of the quadrupole plus toroidal components thus provides another channel for the anti-glitch 
to proceed, one which requires less initial toroidal field strength and less reduction in toroidal field strength, compared to the 
decay of only the toroidal component in a purely dipolar configuration. 

It may be instructive to ex plore some other possible re transitions. Because the volume of the torus changes with re 
(Mastrano, Suvorov. & Melatos 20151) . we find it easier to express our configurations in terms of A, the ratio of the poloidal 
field energy to total field energy. In Fig. 4, we plot some combinations of re/ and A/ (final A) which yield Av/v = —6.3 x 10 - ' 
via equation (10), given initial values A; = 0.1 and re; = 1 (triangles), 0.8 (pluses), and 0.6 (crosses). There is a clear trend that 
lower re/ (for each given re;) and higher re; (for each given re/) need higher A/ (i.e., a greater decrease of toroidal component). 
Figure 4 therefore shows more clearly how the poloidal quadrupole component tends to deform the star into a more oblate 
shape. 


3 Strong higher-order multipoles in magnetars, much stronger than the dipole component (assumed to be the component most responsible 
for spin down), are not ruled out completely. Analysis of the X-ray spectrum of SGR 041 8+5729 suggests that the surface field strength 
is ~ 10 2 stronger than the spin-down-inferred field strength (IRea et al.ll20 1 (ll : iGii ver. Gogtis. fc Ozell l2011). 
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Figure 3. Initial and final toroidal field combination s consistent with the observed an ti-glitch in the AXP IE 2259+586, i.e., solutions 
to equation (10) for /±.v/v = —6.3 x 10 — 7 [cf. Fig. 1 of lGarcfa Sz Ranea-Sandovall d2015h l as a function of initial volume-averaged toroidal 
field strength (-BJ), for Bq = 2.95 x 10 13 G, M* = 1.4Mq, and i?* = 10 4 m. The star has a quadrupole component, which decreases in 
magnitude post-anti-glitch. We keep a = —5 and re.,; = 0.8 (initial re). The three curves represent different final re values: rey = 0.6 (solid 
curve), rey = 0.4 (dashed curve), and rey = 0.1 (dotted curve). The thin solid curve is the identity function, for reference. 


Table 1. Values of the parameters a, b, c, and d in the ellipticity and energy formulae, equations (14) and (16). 
The second and third columns give the values of a and b for different values of re, given in the first column. The 
fourth and fifth columns give the values of c and d for transitions from re = 0.8 to the values of re given in the 
first column. 


Av 

00 

1 

o 

<3 

b 

c 

d 

0.1 

9.37 x 10’ 1 

0.394 

4.90 

0.63 

0.4 

9.96 x 10- 1 

0.236 

4.24 

0.73 

0.6 

1.19 

0.149 

3.64 

0.85 

0.8 

1.77 

0.080 

- 

- 


4 ENERGETICS 


In this section, we calculate the magnetic energy difference between the pre-anti-glitch and post-anti-gli tch states and c ompar e 
it to the observed energy output during the 2012 April 21 burst, which the anti-glitch accompanied ( Archibald et al. 201+ 
For the dipolar case discussed in Section 2, the change in magnetic energy A E m can be written as a function of stellar radius 
and {Bt) as 


A Em = 6.59 x 10 4 


( R, 


V10 4 m 


(B 


f\ 


10 15 G 


For the dipole-plus-quadrupole case described in Section 3, we find 


(Bl) 

10 15 G 


erg. 


(15) 


A E m = c x 10 


f R* 


Vio 4 


AMI )\ d (A\ 

io i5 g ) Go i 5 g; 


erg, 


(16) 


with the dipole surface polar field strength kept constant at 2.95 x 10 13 G. The dimensionless factors c and d are given in 
Table 1 for three representative transitions from Ki = 0.8 to rey = 0.1, 0.4, and 0.6. 

During the anti-glitch, the Fermi and Swift satellites detected an energy release of E ~ 10 38 erg in the 10-10 3 keV band 


and E ~ 10 41 erg in the 2-10 keV band i Folet, et al .11201 j : Archibald et ahlhoTlI : Huang fc Genel2014l : Garcia fc Ranea- Sandoval 


20151) . We obtain, for bot h the dipole and dipole-plus-quadrup ole configurations, A E m ~ 10 4t> —10 4r erg, similar to the most 
energetic SGR giant flare (Palmer et al. 20051 : Mereghetti 20081 ). If the anti-glitch was indeed due to a sudden e change, this 
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Figure 4. Combinations of Kf and A j (final values of k and A) which produce the anti-glitch magnitude consistent with that observed 
in the AXP IE 2259+586 = —6.3 X 10~ 7 ) for a fixed A; (initial A). The star is in a dipole-quadrupole twisted torus configuration, 

with Aj = 0.1 (this initial value of A is indicated by the horizontal dashed line, for reference), and Kj = 1 (triangles), 0.8 (pluses), and 0.6 
(crosses). The dipole surface polar field strength is kept constant at 5.9 X 10 13 G (inferred from spin down) and we assume a canonical 
neutron star mass (1.4 Mq) and radius (10 4 m). 


suggests two possibilities: (1) that, as proposed by iGarcfa fe Ranea-Sandova] ( 2015 ). while the crust readjustment process 


(observed as the anti-glitch) is instantaneous [< 1 d ( [Archibald et al] 120131 )]. the magnetic field decay that leads to it is 


gradual, with time-scale ~ 10 5 yr; or (2) that the magnetic energy is released instantaneously, but is inefficiently (~ 10~ 5 ) 
converted into observed electromagnetic energy. As currently there has only been one confirmed anti-glitch, we cannot draw 
a conclusion yet. 


5 ANTI-GLITCHES IN OTHER MAGNETARS? 


According to equation (9), we need B s > 10 13 G t o get an anti-glitch of the size observed in AXP IE 2259+586 with Aj = 10 -3 
[the lower bound f or stab ility (Akgunetaf 201 3)] or with = —4 x 10 -6 [the maximum ellipticity that the crust can tolerate 
before breaking ( Horowitz fc Kadaul 200gF ]. Why, then, have we not observed an anti-glitch from other magnetars, many of 
which have higher B p than AXP IE 2259+586? This is a question that requires a more detailed study and a more sophisticated 
model than the ones we discuss in this paper. The answer may simply be that we have not yet observed enough objects for 
long enough, or there may be a more microphysical explanation. Either way, using equations (9) and (10), we can set upper 
limits on the possible sizes of anti-glitches from other magnetars. Assuming that (1) the magnetic field is in a purely dipolar 
twisted torus configuration, (2) that the change in e is caused entirely by a change of A (or, equivalently, ( B t )), like in Section 
2, (3) that £i = —4 x 10 -6 (the maximum e that a neutron star’s crust can support before cracking), (4) that ef = e max , and 
(5) M* = 1.4 Mq and R * = 10 4 m, we write down the maximum possible anti-glitch (A v/v) ma , x and the attendant magnetic 
energy change (AE m ) max for four magnetars which have exhibited glitches (even though glitches and anti-glitches may be 
due to entirely different physical processes) and display them in Table 2. 

As evident from Table 2, the maximum anti-glitches should be detectable easily by X-ray timing experiments similar to 
the ones targeting AXP IE 2259+586. The stronger B p is, the larger the size of the maximum anti-glitch and the associated 
magnetic energy change. 

The maximum anti-glitches in Table 2 are accompanied by large magnetic energy releases, stronger the one detected d ur- 
ing AXP IE 22 59+586 anti-glitch, stronger than the 2004 December 27 giant flare of SGR 1806—20 (~ 10 46 erg) ( Palmer et al. 


20051 : Mereghettil 200^) . If we assume the magnetic field change is instantaneous and that ~ 10“° of the released magnetic 
energy is converted into radiation as for AXP IE 2259+586 (as discussed in Section 4), then we can expect the maximal 
anti-glitches from these sample magnetars to be accompanied by giant flares with energies ~ 10 44 -10 4S erg. If these puta¬ 
tive large anti-glitches are not accompanied by energy releases of these magnitudes, we can conclude that, as proposed by 
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Table 2. Maximum anti-glitch magnitudes (Ai'/z') max and the associated magnetic energy change (AE m ) max 
for four glitching magnetars (B p taken from the McGill magnetar catalog), predicted by assuming that the 
magnetic field is in a purely dipolar twisted torus configuration, that the change in e is caused entirely by a 
change of A, like in Section 2, and that one has e = —4 X 10“ 6 initially [the maximum ellipticity that the crust 
can support before breaking (jllorowit.z &; Kadaull200~9l )l. that £f = 0, and M* = 1.4 Mq and f?* = 10 4 m. 


Name 

B p (10 14 G) 

(A(z/iz) max (lO- 6 ) 

(AE m ) max (10 48 erg) 

4U 0142+61 

1.3 

-2.7 

-9.2 x 10 _1 

IE 1048.1-5937 

3.9 

-3.0 

-1.0 

1RXS J170849.0—40091061 

4.7 

-3.2 

-1.1 

IE 1841-045 

7.0 

-3.9 

-1.3 


Garcia fe Ranea-Sandoval ( 20151) . the magnetic fields must decay gradually before the crust cracks and readjusts, or that the 


electromagnetic conversion efficiency varies substantially across the population. 


6 DISCUSSION 

In this paper, we refine the model proposed by Garcia fe Ranea-Sandoval ( 20151 ) to explain the anti-glitch of the AXP IE 
2259+586 in terms of a sudden change in ellipticity e. We also show how the observed anti-glitch can also occur via an 
alternative magnetic channel, namely the decay of a higher-order multipolar field (in addition to a decay of the toroidal 
component). We calculate the change in e self-consistently for a dipole twisted torus in Section 2 and a dipole-quadrupole 
twisted torus in Section 3. We show that, by neglecting the contribution of the poloidal component to the Lorentz force, one 
overestimates the prolateness of the star (for example, for 1.9 x 10 ls G < (Bt) < 4.9 x 10 15 G and B p = 5.9 x 10 13 G, we 
find e/eGRS ~ 0.15), which leads to an underestimate of the required (B\) and the required change in field strength (Figs. 1 
and 2). For both dipole and dipole-quadrupole twiste d tor i, we find that onl y objects w ith B p > 10 13 G (e.g., magnetars) can 
match the observed /S.v/v, confirming a conclusion of Garcia fc Ranea-S ando val (|2015|). This means that we are not likely to 
observe a similar anti-glitch in a radio pulsar and that, despite superficial similarities, a magnetar anti-glitch may have an 
entirely different physical origin from a magnetar or pulsar glitch. 

If the interior of a neutron star consists of superconducting protons (lLanderl2013l ) or quarks dGlampedakis. Jones, fc Samuelssonl 
20121). t hen [e]_ is raised by a factor of ~ H c i/(B ), where H c 1 ~ 10 16 G is the lower superconductivity critical field 


I G lampedakis . A ndersso n. fe Sa muelss on 120111 ). In AXP IE 2259+586, for example, a superconducting interior can raise 
|e| by a factor of ~ 10, lowering the required (Bt) and A (Bt) by a factor of ~ 0.3. Future anti-glitches, if observed in a star 
with a lower dipole field inferred from spin down, may be taken as evidence of superconducting interiors. 

For e to change by the amount needed to give the observed tSv/v in AXP IE 2259+586, the magnetic energy change 
is larger by 5 or 6 orders of magnitude than the observed outburst energy. In Section 5, we calculate the upper limits of 
anti-glitch magnitudes and energies in four other glitching magnetars (assuming for definiteness that they start at the edge of 
stability, A = 10 -3 ). We find (Aiz/j/) max ~ 10 -J -10 -3 , higher than the observed AXP IE 2259+586 anti- glitch, and m agnetic 
ener g y releases o f order 10 49 -10 so erg, higher than the most energetic magnetar giant flare of ~ 10 46 erg I Palmer et all 2005) : 


Mereghetti 2008 ). If these anti-glitches are observed in the future, we can use the energy of the accompanying burst/flare (if 


any) to help conclude whet her the field reconfiguration is instantaneous but with an inefficient radiative conversion or the 
field decay is a slow process I Garcia fe Ranea-Sandoval 20151 ). Lander et al. 1 20151 ) recently set an upper limit of 4 x 10 46 erg 
on energy released by crustal fracture during a magnetar flare (largely independent of magnetic field strength), assuming the 
fracture extends to the base of the crust. Therefore, a high-energy flare can, in principle, be powered solely by a ‘crustquake’. 
Differentiating the energy contributions from crust and field may be difficult in practice, but if a flare is accompanied by an 
anti-glitch, a significant fraction of the observed energy may be due to field reconfiguration. 

The anti-glitch in AX P IE 2 259+586 has also been interpreted as evidence that AXPs and other magnetars are surrounded 
by fallback matter ( Katz 20131 ). Huang fe Gene ( 2014 ) p roposed that the anti-glitc h is due to a collision between AXP IE 
2259+586 and a small solid body with a mass ~ 10 18 kg. IKanto r fc Gusakov ( 20141 ) suggested instead that the anti-glitch is 
a consequence of neutron star superfluidity, just like a normal pulsar glitch, due to a velocity lag between the superfluid and 
the crust. We do not consider these alternative scenarios in this paper. 

We also do not model the trigger and time-scale of the anti-glitch. Like [Garcia fe Ranea-Sandoval (2015), we implicitly 
assume that the anti-glitch (and the outburst it accompanied) occurred, when the magnetar crust cra cked d ue to some internal 
instability, built up as the internal toroidal magnetic field decayed (but see also Sections 4 and 5). Link |2014), in contrast, 
concluded that magnetar bursts and flares must be ca used by a r elaxat ion of the external field, to be consistent with the rise 
times of the accompanying quasi-periodic oscillations. Lvutikov ( 2015) ) similarly concluded that a magnetar burst/flare can 


© ? RAS, MNRAS 000, [THU] 












































































AXP IE 2259+586 antiglitch and its internal magnetization 11 


proceed without any internal instability leading to crust cracking. A simulation of the field decay process, the resultant crust 
cracking, and a calculation of the time-scales involved are reserved for future work. 
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